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ABSTRACT 

In this paper we analyse the effects of an axial electric field on 
helical flow of incompressible cholesteric liquid crystals between two 
rotating coaxial circular cylinders having different axial velocities. The 
system of governing differential equations for the flow have been obtained 
and is then solved numerically by the finite difference method. The results 
obtained are then graphically represented. It has been observed that the 
velocities and the orientations of the molecules decreases with the 
increase in the intensity of electric field. 
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1. INTRODUCTION 

The past decades have witnessed the 
modeling of static versions of liquid crystals 
by Oseen 15 and Zocher 22 . This work was 
subsequently modified by Frank 6 by taking into 
account of the physical configurations of liquid 
crystals. Extending the work Ericksen 4 ' 5 refor- 
mulated this theory to describe the dynamical 
and hydrostatic behavior of liquid crystals. On 
the basis of this work Leslie 1011 proposed a set 
of constitutive equations for studying the 



dynamic behavior of cholesteric liquid crystals 
in the presence of thermal effects. Later on Currie 3 
discussed the helical flow of nematic liquid 
crystals. In order to examine the effects of electric 
field on liquid crystals, Skarp et al. 19 have proposed 
the modified expressions for the electromag- 
netic part of stress tensor and Helmholtz free 
energy function. 

Normally liquid crystals are classified 
into three categories : nematic, cholesteric and 
smectic. Amongst these the nematic liquid 
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crystals possess a high degree of long range 
orientational order of the molecules, but not 
comprising of the long range translational order, 
while the cholesteric mesophase resembles the 
nematic type of liquid crystal except with a 
difference that it is composed of optically active 
molecules. Consequently the theory of cholesteric 
liquid crystals 11 is much more complex than the 
theory of nematic liquid crystals 12 . Due to this 
less attention has been paid to the investigations 
concerning cholesteric liquid crystals. 

The pitch of cholesterics is very sensitive 
to the external fields but is important in the 
applications to the display devices. Also since 
Liquid crystal films are inexpensive can work under 
the low voltage and low power and therefore 
interest in cholesterics has been rapidly 
growing. In the last few years Morais et a/. 13 , 
Prishchepa et al. 17 , Rodriguez-Rosales et a/. 18 
and Nakata et al. 14 have found different types 
of applications of liquid crystals. Also there has 
been numerous experimental and theoretical 
investigations 1 ' 2 ' 7 ' 9 ' 20,21 concerning cholesteric 
liquid crystals. 

In this paper, we propose to apply a 
numerical approach to evaluate the effects of 
electric field on the helical flow of incompressible 
cholesteric liquid crystals between two coaxial 
circular cylinders of infinite extent. Here both 
the cylinders are rotating with different angular 
velocities and are simultaneously moving along 
their common axis with different axial 
velocities. The boundary conditions used in this 
problem are adopted to be the same as used by 
Paria et al} 6 elsewhere. The electric field has 
been directed along the common axis of 
cylinders. 



2. Basic Equations : 

The basic equations for flow of incom- 
pressible cholesteric liquid crystals as used by 
Leslie 11 are 

(i) the equation of continuity : 

n i = o (2.i) 

(ii) the equation of motion of incompressible 
liquid crystals with respect to mechanical 

velocity y'- 

Dv, 

P-^~ =P p i + °ju (2-2) 

(iii) the equation of motion of director d to 
describe the motion of liquid crystal molecules : 

D ^ cl 

Pi— T~ = PiG, + gt+xju. (2.3) 
Dt 

Here the symbols like y , p, F t , a ji, pi, 
D 

d , Gi, gi, Tiji, and are used to denote the 

velocity vector, the uniform density, the body 
force per unit mass, the stress tensor, an inertial 
coefficient , the director , the extrinsic director 
body force per unit mass, the intrinsic director 
body force per unit volume , the director stress 
tensor and the material time derivative 
respectively. 

Following Skarp et al} 9 the above 
basic equations in the presence of electric fields 
are supplemented by the following set of 
constitutive equations for the stress tensor ay,-, 
the director stress tensor and the intrinsic 
director body force g ; 

Oji = -pdij + ffjt, (2.4) 
Xji = Pjdi+ p-^j— + cce ijk d k , (25) 
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g, = ;•'/, /'';'/,/ P~^- -oce ijk d k ,j + f,.(2. 6 ) 

Here the symbols p, F, a, ft ,J, ^ji and 

gi represent the pressure, the Helmholtz free 
energy per unit mass, a material coefficient, an 
arbitrary vector, the director tension, the distortion 
stress and the non-equilibrium part of extra 
intrinsic director body force respectively The 
expression for the distortion stress and the 
contributions to non-equilibrium part of extra 
intrinsic director body force in the absence of 
thermal effects are 



a H = - P 



DF 



dd 



d ki+ ae jk P ( d pdi),k 



h = KN i + k 2 A ik d k , 



(2.7) 
(2.8) 



The electromagnetic part of stress tensor 
<J ^ em in terms of electric displacement A and 

electric field component E t defined by Skarp 
et al. 19 is 



_ em 



DiEj+DjEi-faE^q 



(2.11) 



where 

D i ={£ ± S ij +£ a d i d j )E j . (2.12) 

Here s a and e ± denote the dielectric anisotro- 
pics of the molecules. 

The modified form of Helmholtz free 
energy Fused by Skarp et al. 19 for cholesteric 
liquid crystals is 

IpF = a l (d ij ) + a 1 (x + d l e ljk d kj ) + a^d { d j d t , j d k j + 
(« 2 + « 4 Mi, j dj, i ~ (d,_ if]- £>,£'„ (2. 1 3) 



where the symbols <T j ; and <x J( em denote 

the non-equilibrium part of stress tensor and the 
electromagnetic part of stress tensor respectively. 
Further the non- equilibrium parts of stress 
tensor is given by 

Sji = mdkdpAfydidj + m N i d j + H3 N j d i 
+ n 4 A y + p 5 A ik d k d j 
+ V 6 A jk d k d i . ( 2 -9) 

Here the material coefficients fii and 2, are 
related to each other by the relations 

K =^2-^3' X 2 =M5~M 6 and 

2A U = Vij + vjj, - 



2W<j = Vij - v p , 



(2.10) 



>- 



where a , and t are the material coefficients. 

All the material coefficients occurring in the 
above equations are taken as constant. 

3. Statement of Problem : 

In this paper we wish to examine the 
helical flow of an incompressible cholesteric 
liquid crystal with director of unit magnitude 
between two infinite coaxial circular cylinders 
in the presence of a electric field. The cylinders 
are rotating and also moving with different axial 

velocities and the electric field E is acting along 
the common axis of the cylinders. 

We choose a system of cylindrical polar 

co-ordinates {r,9,z) such that the z-axis 
coincides with the common axis of the cylinders. 



The outer and inner cylinders of radii r 2 and rj, 
(r 2 > n) are rotating with angular velocities 002 
and CO] and are moving with axial velocities u 2 
and ui respectively. 

We shall analyse the solution of 
equations (2.1) to (2.3) in the form [16] 

d r = sin <j>{r\do = cos ^(r)sin y/(r\d z = cos ^(r)cos y/{r), 
v r = 0,v e = rco(r), v z = u(r) (3.1) 

where the angle <f> is being used for 
the translational motion of the molecules while 
the twist angle (//represents the spinning motion 
of the molecules of liquid crystal. 

Following Paria et al. 16 the boundary 
conditions for the velocities and the orientation 
of molecules are taken as 



u = u x , 



02=0, 



co = co l 



a = ai 0 



¥\ = ¥o 
¥2 = ¥0 



at r = r x , 
at r = r 2 , 

at r = r x , 
at r = Tt. 



(3.2) 



(3.3) 



4. Formulation of differential equations: 

We observe that the equation of 
continuity (2.1) is satisfied by our choice of 
velocity components. The equations of motion 
(2.2) and (2.3) in our case take the forms 



A rn rr) 

dr 
dr 
dr 



+ 7tn 



+ s z =o 



+ g r +p l co sin^ = 0 > (4.4) 
+ g e + Pico 1 cos (j) sin 1// =0 (45^ 
(4.6) 



Equation (4.1) with the help of equation 
(2.4) gives 



P = P 0 + 



2 da rr a rr - d g 

pra + rjL - + — 

dr r 



dr 



(4.7) 



where P 0 is the pressure on the inner cylinder 
and 



a r =P + a a , <7„„=P + <7 ( 



'00 



Also equations (4.2) and (4.3) with the 
help of the equations (2.4), (2.7) , (2.9) , (2. 10), 
(2.11) and (2.12) on integration leads to the 
following equations in t, and 77 as 



( . . . ) --[ ( 2 ( )si 



sin cos } 



dr 
d(ra r9 ) 



dr 

d(r<J rz ) 
dr 



+ a 6r =0 



(4.1) 
(4.2) 
(4.3) 
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^ h 
-2fl^ + 



2/z 



2a 4 (am + ariyRj 



(e 4 +e 3 +e 6 ) 
2d(am + an)(e 3 + e 6 ) 2Elg 4 (am + an)R i 2 



e 4 [e 4 +e 3 +e 6 



) 



2d(am + an) 2 , n r< i\ 

+ — V,- -4£l£ 4 fl,. 0 f = 0, (5.7) 



a 2 \R{ 



Vi+h ~ 2 Vi +Vi-h I , # ( Vi+fc 



2h 



-a 2 yr,. ^-2£lg 4J R,.V ; =0.(5.8) 

Here a 3 =e 4 +e 3 +e 6 , a 4 =e 3 +e 6 , n=e 4 , 

11 k k 

8\ = = ,82 = 



2 



2m 1(Ui 



2a, 



'1 



Since the orientations 0 ; and ^, are involved 

Uj / 

in the equations of the velocities / u and 



CO; 



, we shall now determine <f> i and Wi first 



by using the boundary conditions (3.3) in the 
equations (5.7) and (5.8) by making use of the 
Cramer's rule. Thereafter with the help of these 
values and the boundary conditions (3.2) the 

axial velocities and the angular velocities 



CO; 



'co x 



are computed with the help of equations 



(5.5) and (5.6). Throughout the above calcula- 
tions it has been assumed that 



%= 2 '%r 2 < v/o=1 ' h=0 - 2 ' r > =i > 

where j = 1, i and/ = 1.2, 1.4, 1.6, 1.8. 

From these results it is observed that the electric 
field parameters El occur in the equations of 

Uj, ax, <f>i and (//,-. 

6. RESULT AND DISCUSSION 

We shall now numerically compute the 
velocities and the orientations by using the 
equations (5.5), (5.6), (5.7) and (5.8). These 
numerical values are then plotted in figs. (1) to 
(4) versus the radial distance (r/ri) for different 
values of the electric field parameter El , where 
r/ri is varied from 1 to 2 with a step size of 0.2. 

From figs. (1) and (2) which illustrate 
the axial (u/ui) and angular velocities (co/coi) 
versus r/ri profiles for different values of the 
electric field parameter El, it is observed that 
the increase of radial distance from the inner 
to the outer cylinder leads to sinusoidal form 
of behaviour of u/ui and co/coi upto the first half 
and subsequently it repeats again the sinusoidal 
form of behaviour with a little rise in the first 
contour whereas a sharp decrease is indicated in 
the second contour. However both type of 
velocities decrease continuously with the increase 
in the magnitude of El. 

Fig. (3) and fig. (4) illustrate the effects 
of the electric field parameter El on the orientation 
profiles for <f> and y/ respectively. It is being 
noticed that the profiles of <f> and y/ decrease 
continuously with the increase in the values of 
El. However for a fixed value of El the 
orientation <f> increases upto mid of the distance 
and then it decreases, but in case of orientation 
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y/ the graph decreases first upto the mid point 
and afterwards it rises upto the other end. 

7. CONCLUSION 

The flow of cholesteric liquid crystals 
between two co-axial circular cylinders under 
external influence is of engineering and industrial 
importance. Also the change in orientation 
caused by electric field effects in cholesteric liquid 
crystals can be used in liquid crystals display 
devices. Scientists and engineers are using 
cholesteric liquid crystals in a variety of applica- 
tions because external perturbations cause 
significant changes in the macroscopic properties 
of the liquid crystals. The electric field is capable 
of inducing these changes. However the 
magnitude of the field as well as the speed at 
which the molecules align are also the important 
characteristics that industry has to deal with. 
We, therefore, hope that the present investigations 
possess the possibility of such applications. 
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